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1. Introduction 



In this chapter, we will be concerned with hadronic interactions in the nuclear medium, 
particularly imder conditions of extreme densities such as those encountered in some compact 
astrophysical objects. This issue goes to the very core of nuclear physics. In fact, our present 
knowledge of the nuclear force in free space is, in itself, the result of decades of struggle 
jMachleidt |1989| which will not be reviewed in this chapter. The nature of the nuclear force 



in the medium is of course an even more complex problem, as it involves aspects of the force 
that cannot be constrained through free-space nucleon-nucleon (NN) scattering. Predictions 
of properties of nuclei are the ultimate test for many-body theories. 

Nuclear matter is a convenient theoretical laboratory for many-body theories. By "nuclear 
matter" we mean an infinite system of nucleons acted on by their mutual strong forces and no 
electromagnetic interactions. Nuclear matter is characterized by its energy per particle as a 
function of density and other thermodynamic quantities, as appropriate (e.g. temperature). 
Such relation is known as the nuclear matter equation of state (EoS). The translational 
invariance of the system facilitates theoretical calculations. At the same time, adopting what 
is known as the "local density approximation", one can use the EoS to obtain information on 
finite systems. This procedure is applied, for instance, in Thomas-Fermi calculations within 
the liquid drop model, where an appropriate energy functional is written in terms of the EoS 
jFurnstahH[2002}|Oyamatsu et al.|[l998,:,Sammarruca & Liu„2009, |. 

Isospin-asymmetric nuclear matter (lANM) simulates the interior of a nucleus with unequal 
densities of protons and neutrons. The equation of state of (cold) lANM is then a function of 
density as well as the relative concentrations of protons and neutrons. 



The recent and fast-growing interest in lANM stems from its close connection to the physics 
of neutron-rich nuclei, or, more generally, isospin-asymmetric nuclei, including the very 
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"exotic" ones known as "halo" nuclei. At this time, the boundaries of the nuclear chart are 
uncertain, with several thousands nuclei believed to exist but not yet observed in terrestrial 
laboratories. The Facility for Rare Isotope Beams (FRIB) has recently been approved for design 
and construction at Michigan State University. The facility will deliver intense beams of 
rare isotopes, the study of which can provide crucial information on short-lived elements 
normally not found on earth. Thus, this new experimental program will have widespread 
impact, ranging from the origin of elements to the evolution of the cosmos. It is estimated 
that the design and construction of FRIB will take ten years. In the meantime, systematic 
investigations to determine the properties of asymmetiic nuclear matter are proliferating at 
existing facilities. 

The equation of state of lANM is also the crucial input for the structure equations of compact 
stars, and thus establishes the connection between nuclear physics and compact astrophysical 
systems. It is the focal point of this chapter to present and discuss our approach to the 
development of the EoS of nuclear and neutron-rich matter, with particular emphasis on the 
latter and its relation to the properties of neution stars. 

The chapter will articulate through the following sections: In Section 2, we present a 
brief review of facts and phenomenology about lANM. We then proceed to describe our 
microscopic approach to calculate the energy per particle in lANM (Section 3) and show the 
corresponding predictions. Section 4 will be dedicated to a review of neutron star structure 
and available empirical constraints. Microscopic predictions of neutron star properties will be 
presented and discussed there. Section 5 contains a brief discussion on the topic of polarized 
lANM. The chapter will end with our conclusions and an outlook into the future (Section 6). 

2. Facts about isospin-asymmetric nuclear matter 

Asymmetiic nuclear matter can be characterized by the neution density, pn, and the proton 
density, pp, defined as the number of neutions or protons per unit of volume. In infinite 
matter, they are obtained by summing the neution or proton states per volume (up to their 
respective Fermi momenta, kp or kp) and applying the appropriate degeneracy factor. The 
result is 

with i = n or p. 

It may be more convenient to refer to the total density p = pn + pp and the asymmetry (or 
neutron excess) parameter a = £2-_ft Clearly, a=0 corresponds to symmetiic matter and oc=l 
to neutron matter. In terms of a and the average Fermi momentum, kp, related to the total 
density in the usual way, 

the neutron and proton Fermi momenta can be expressed as 

fc? = A:f(l + a)i/3 (3) 

and 

fcP=fcf(l-«)i/3, (4) 
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respectively. 

Expanding the energy per particle in lANM with respect to the asymmetry parameter yields 

where the first term is the energy per particle in symmetric matter and the coefficient of the 
quadratic term is identified with the symmetry energy Csym ■ In the Bethe-Weizsacker formula 
for the nuclear binding energy it represents the amount of binding a nucleus has to lose when 
the numbers of protons and neutrons are unequal. A typical value for egym at nuclear matter 
density (po) is 30 MeV, with theoretical predictions spreading approximately between 26 and 
35 MeV. 

To a very good degree of approximation, the energy per particle in lANM can be written as 

e{p,a.) K. eQ{p)+esym{p)a^. (6) 

The effect of a term of fourth order in the asymmetry parameter (C (a^)) on the bulk properties 
of neutron stars is very small, although it may impact the proton fraction at high density. 

Equation |6| displays a convenient separation between the symmetric and aymmetric parts 
of the EoS, which facilitates the identification of observables that may be sensitive, for 
instance, mainly to the symmetry energy. At this time, groups from GSI jSfienti et al.||2009 



Trautmann et al."20 09|, M SU (fTsanget al.; 2009'l, Italy fGreco'.'2010'l, France, fBorderie & Rivet 



2008), China I F eng| |2010] [Vbn g 2010), and Japan (Isobe, 20111 are investigating the density 



dependence of the symmetry energy through heavy-ion collisions. Typically, constraints are 
extracted from heavy-ion collision simulations based on transport models. Isospin diffusion 
and the ratio of neutron and proton spectra are among the observables used in these analyses. 

These investigations appear to agree reasonably well on the following parametrization of the 
symmetry energy: 



{p) = 12.5 MeV (-^Y^^ + 17.5Mev(-^y', (7) 
\pQ/ \pQ/ 



where po is the saturation density. The first term is the kinetic contribution and 7, (the 
exponent appearing in the potential energy part) is found to be between 0.4 and 1.0. Recent 
measurements of elliptic flows in ^'''Au + ^^''Au reactions at GSI at 400-800 MeV per nucleon 
favor a potential energy term with 7,- equal to 0.9 ± 0.4. Giant dipole resonance excitation in 
fusion reactions jSimenel et al.[[2007) is also sensitive to the symmetry energy, since the latter 
is responsible for isospin equilibration in isospin-asymmetric collisions. 

Isospin-sensitive observables can also be identified among the properties of normal nuclei. 
The neutron skin of neutron-rich nuclei is a powerful isovector observable, being sensitive to 
the slope of the symmetry energy, which determines to which extent neutrons will tend to 
spread outwards to form the skin. 



Parity-violating electron scattering experiments are now a realistic option to determine 
neutron distributions with imprecedented accuracy. The neutron radius of ^'^^Pb is expected 
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to be measured with a precision of 3% thanks to the electroweak program at the Jefferson 
Laboratory, the PREX experiment in particular, just recently completed at Jefferson Lab. This 
level of accuracy could not be achieved with hadronic scattering. Parity-violating electron 
scattering at low momentum transfer is especially suitable to probe neutron densities, as the 
Z" boson couples primarily to neutrons. With the success of this program, reliable empirical 
information on neutron skins will be able to provide, in turn, much needed independent 
constraint on the density dependence of the symmetry energy. 

A measure of the density dependence of the symmetry energy is the symmetry pressure, 
defined as 

where we have used Eq. |6| with Thus, L is sensitive to the gradient of the energy per 
particle in neutron matter (en.m.)- As to be expected on physical grounds, the neutron skin, 
given by 

S=^<rl>- ^< rj > , (9) 
is highly sensitive to the same pressure gradient. 

Values of L are reported to range from -50 to 100 MeV as seen, for instance, through the 
numerous parametrizations of Skyrme interactions, all chosen to fit the binding energies 
and the charge radii of a large number of nuclei ( jLi & Chen 2005^ and references therein). 



Heavy-ion data impose boundaries for L at 85 ± 25 MeV, with more stringent constraints being 
presently extracted. At this time constraints appear to favor lower values of the symmetry 
pressure. In fact, a range of L values given by 52.7 ± 22.5 MeV has emerged from recent 
analyses of global optical potentials ( |Xu et al.[|2010| . 

Typically, parametrizations like the one given in Eq. |7| are valid at or below the saturation 
density, pQ. Efforts to constrain the behavior of the symmetry energy at higher densities 
are presently being pursued through observables such as tt" -/n+ ratio, K+/K° ratio, 
neutron/proton differential transverse flow, or nucleon elliptic flow jKo et al. 2010^ . 



Another important quantity which emerges from studies of lANM is the symmetry potential. 
Its definition stems from the observation that the single-particle potentials experienced by 
the proton and the neutron in lANM, iJn/p/ different from each other and satisfy the 
approximate relation 

U„/p(k,p,a) !i„/p (fc, p,K = 0)± Usym (k, p) a , (10) 
where the +(-) sign refers to neutrons (protons), and 

Usym = . (11) 

(Later in the chapter we will verify the approximate linear behavior with respect to ol 
displayed in Eq. (TO).) Thus, one can expect isospin splitting of the single-particle potential 
to be effective in separating the collision dynamics of neutrons and protons. Furthermore, 
U-symr being proportional to the gradient between the single-neutron and the single-proton 
potentials, should be comparable with the Lane potential jLane 1962| , namely the isovector 
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part of the nuclear optical potential. Optical potential analyses can then help constrain this 
quantity and, in turn, the symmetry energy. 

Because of the fundamental importance of the symmetry energy in many systems and 
phenomena, it is of interest to identify the main contributions to its density dependence. In 
a recent work | |Sammarruca| [201 1 1 we discussed the contribution of the isovector mesons (tt, 
p, and S) to the symmetry energy and demonstrated the chief role of the pion. Note that 
the isovector mesons carry the isospin dependence by contributing differently in different 
partial waves, and that isospin dependence controls the physics of lANM. Hence, we stress 
the relevance of a microscopic model that contains all important couplings of mesons with 
nucleons. 



3. Our microscopic approach to isospin-asymmetric nuclear matter 



3.1 The two-body potential 



Our approach is ab initio in that the starting point of the many-body calculation is a realistic 
NN interaction which is then applied in the nuclear medium without any additional free 
parameters. Thus the first question to be confronted concerns the choice of the "best" 
NN interaction. After the development of Quantum Chromodynamics (QCD) and the 
understanding of its symmetries, chiral effective theories jWeinberg] |1990 ' were developed 



as a way to respect the symmetries of QCD while keeping the degrees of freedom (nucleons 
and pions) typical of low-energy nuclear physics. However, chiral perturbation theory (ChPT) 
has definite limitations as far as the range of allowed momenta is concerned. For the purpose 
of applications in dense matter, where higher and higher momenta become involved with 
increasing Fermi momentum, NN potentials based on ChPT are unsuitable. 

Relativistic meson theory is an appropriate framework to deal with the high momenta 
encountered in dense matter. In particular, the one-boson-exchange (QBE) model has proven 
very successful in describing NN data in free space and has a good theoretical foundation. 
Among the many available OBE potentials, some being part of the "high-precision generation" 
jMachleidt]|2001 Stocks et al. |1994[|Wiringa et aT 19951, we seek a momentum-space potential 



developed within a relativistic scattering equation, such as the one obtained through the 
Thompson (Thompson 19701 three-dimensional reduction of the Bethe-Salpeter equation 
jjSalpeter & Bether^ ,1951j . Furthermore, we require a potential that uses the pseudovector 
coupling for the interaction of nucleons with pseudoscalar mesons. With these constraints in 



mind, as well as the requirement of a good description of the NN data, Bonn B (Machleidt 
1989| is a reasonable choice. As is well known, the NN potential model dependence of 



nuclear matter predictions is not negligible. The saturation points obtained with different 
NN potentials move along the famous "Coester band" depending on the strength of the tensor 
force, with the weakest tensor force yielding the largest attraction. This can be imderstood in 
terms of medium effects (particularly Pauli blocking) reducing the (attractive) second-order 
term in the expansion of the reaction matrix. A large second-order term will undergo a large 
reduction in the medium. Therefore, noticing that the second-order term is dominated by the 
tensor component of the force, nuclear potentials with a strong tensor component will yield 
less attraction in the medium. For the same reason (that is, the role of the tensor force in 
nuclear matter), the potential model dependence is strongly reduced in pure (or nearly pure) 
neutron matter, due to the absence of isospin-zero partial waves. 
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Fig. 1. Contribution to the NN interaction from virtual pair excitation. Upward- and downward-pointing 
arrows represent nucleons and antinucleons, respectively. Dashed lines denote mesons. 



Fig. 2. Three-body force due to virtual pair excitation. Conventions as in the previous figure. 



Already when QCD (and its symmetries) were unknown, it was observed that the contribution 
from the nucleon-antinucleon pair diagram, Fig. [T] becomes imreasonably large if the 
pseudoscalar (ps) coupling is used, leading to very large pion-nucleon scattering lengths 
jBrown[ |1979). We recall that the Lagrangian density for pseudoscalar coupling of the nucleon 
field (xp) with the pseudoscalar meson field {(p) is 



(12) 



On the other hand, the same contribution, shown in Fig. [Tj is heavily reduced by the 
pseudovector (pv) coupling (a mechanism which became known as "pair suppression"). The 
reason for the suppression is the presence of the covariant derivative at the pseudovector 
vertex, 

/ps 



pv 



(13) 



which reduces the contribution of the vertex for low momenta and, thus, explains the small 
value of the pion-nucleon scattering length at threshold jBrown 1979^ . Considerations 
based on chiral symmetry ( [Brown 1979^ can further motivate the choice of the pseudovector 
coupling. 

In closing this section, we wish to highlight the most important aspect of the ab initio approach: 
Namely, the only free parameters of the model (the parameters of the NN potential) are 
determined by fitting the free-space NN data and never readjusted in the medium. In other 
words, the model parameters are tightly constrained and the calculation in the medium is 
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parameter free. The preserice of free parameters in the medium would generate effects and 
sensitivities which are hard to control and interfere with the predictive power of the theory. 

3.2 The Dirac-Brueckner-Hartree-Fock approach to symmetric and asymmetric nuclear 
matter 

3.2.1 Formalism 

The main strength of the DBHF approach is its inherent ability to account for important 
three-body forces through its density dependence. In Fig. |2] we show a three-body force 
(TBF) originating from virtual excitation of a nucleon-antinucleon pair, known as "Z-diagram". 
Notice that the observations from the previous section ensure that the corresponding diagram 
at the two-body level. Fig. [T] is moderate in size when the pv coupling, Eq. [13] is used. The 
main feature of the DBHF method turns out to be closely related to the TBF depicted in Fig.|2] 
as we will argue next. In the DBHF approach, one describes the positive energy solutions of 
the Dirac equation in the medium as 

where the nucleon effective mass, m*, is defined as m* = m + Ug, with Ug an attractive 
scalar potential. (This will be derived below.) It can be shown that both the description of a 
single-nucleon via Eq. l |14| and the evaluation of the Z-diagram, Fig. |2j generate a repulsive 
effect on the energy per particle in symmetric nuclear matter which depends on the density 
approximately as 

^E.{^Y\ (15) 

and provides the saturating mechanism missing from conventional Brueckner calculations. 
(Alternatively, explicit TBF are used along with the BHF method in order to achieve a similar 
result.) Brown showed that the bulk of the desired effect can be o btained as a lowes t order 



(in p'^/m) relativistic correction to the single-particle propagation | |Brown et al. 



1987'. With 



the in-medium spinor as in Eq. (14) , the correction to the free-space spinor can be written 
approximately as 





m*(P,A)-m(p,A) w ir-p,j ] Xa, (16) 

where for simplicity the spinor normalization factor has been set equal to 1, in which case 
it is clearly seen that the entire effect originates from the modification of the spinor 's lower 
component. By expanding the single-particle energy to order Ug, Brown showed that the 



correction to the energy consistent with Eq. Jl6| can be written as ( ^jf He then proceeded 



to estimate the correction to the energy per particle and foimd it to be approximately as given 
in Eq. jTs). 

The approximate equivalence of the effective-mass description of Dirac states and the 
contribution from the Z-diagram has a simple intuitive explanation in the observation that 
Eq. |T4) , like any other solution of the Dirac equation, can be written as a superposition of 
positive and negative energy solutions. On the other hand, the "nucleon" in the middle of 
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the Z-diagram, Fig. |2j is precisely a superposition of positive and negative energy states. In 
summary, the DBHF method effectively takes into account a particular class of TBF, which are 
crucial for nuclear matter saturation. 

Having first summarized the main DBHF philosophy, we now proceed to describe the DBHF 
calculation of lANM jAlonso & Sammarruc a 2003 S ammarruca[ |2010| . In the end, this will 
take us back to the crucial point of the DBHF approximation, Eq. jl4^ . 

We start from the Thompso n ^Thompson) |1970| rela tivistic three-dimensional reduction of 



the Bethe-Salpeter equation (Salpeter & Bether 1951^. The Thompson equation is applied to 



nuclear matter in strict analogy to free-space scattering and reads, in the nuclear matter rest 
frame, 

+ ./ J2^^,^^' ^) Ef4 ie*l^t.iP,Kf '^ ^' ^^^^ ' ^''^ 

where gij is the in-medium reaction matrix {ij=nn, pp, or np), and the asterix signifies that 
mediimi effects are applied to those quantities. Thus the NN potential, v*j, is constructed 



in terms of effective Dirac states (in-medium spinors) as explained above. In Eq. 1 17 , q, q' , 



and K are the initial, final, and intermediate relative momenta, and E* = {m*Y + T^^- The 
momenta of the two interacting particles in the nuclear matter rest frame have been expressed 
in terms of their relative momentum and the center-of-mass momentum, P, through 

P = ki+k2 (18) 

and 

(19) 

The energy of the two-particle system is 

e*. (P, K) =e*{P,K)+ e* (P, K) (20) 

and (e,^)o is the starting energy. The single-particle energy e* includes kinetic energy and 
potential energy contributions. The Pauli operator, Q,y, prevents scattering to occupied nn, pp, 
or np states. To eliminate the angular dependence from the kernel of Eq. | (T7| , it is customary 
to replace the exact Pauli operator with its angle-average. Detailed expressions for the Pauli 
operator and the average center-of-mass momentum in the case of two different Fermi seas 
can be found in Ref. jAlonso & Sammarruca[[2003| . 

With the definitions 

m* m* 

and 
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one can rewrite Eq. |T7) as 

G,0,c},P,(£*j)o) = V*jiq\ci) 



■I 



(2^)' ' '{£*Ao-£UP,K) 



,7)0 -e,/ 



which is formally identical to its non-relativistic counterpart. 

The goal is to determine self-consistently the nuclear matter single-particle potential which, 
in lANM, will be different for neutrons and protons. To facilitate the description of the 
procedure, we will use a schematic notation for the neutron/ proton potential. We write, for 
neutrons, 

!i„ = !i„p + Unn , (24) 

and for protons 

!ip = !Jp„ + iJpp , (25) 

where each of the four pieces on the right-hand-side of Eqs. i24j25j signifies an integral 
of the appropriate G-matrix elements (nn, pp, or np) obtained from Eq. p3) . Clearly, the 
two equations above are coupled through the wp component and so they must be solved 
simultaneously. Furthermore, the G-matrix equation and Eqs. j24|25| are coupled through 
the single-particle energy (which includes the single-particle potential, itself defined in terms 
of the G-matrix). So we have a coupled system to be solved self-consistently. 

Before proceeding with the self-consistency, one needs an ansatz for the single-particle 
potential. The latter is suggested by the most general structure of the nucleon self-energy 
operator consistent with all symmetry requirements. That is: 

Uiip) = Us,iip) + 7oU°v,dP) -y-P UvAp) . (26) 

where (Js_, and are an attractive scalar field and a repulsive vector field, respectively, 
with iiy ; the timelike component of the vector field. These fields are in general density and 
momentum dependent. We take 

W,(p)«LIs,,(p) + 7o</(p), (27) 

which amounts to assuming that the spacelike component of the vector field is much smaller 
than both Lig , and Uy ^. Furthermore, neglecting the momentum dependence of the scalar and 
vector fields and inserting Eq. j27| in the Dirac equation for neutrons/ protons propagating in 
nuclear matter, 

(7;,p''-?«,-W,(p))m,(P.A)=0, (28) 
naturally leads to rewriting the Dirac equation in the form 

(7,,(pf')*-'«n"KP.A) = 0, (29) 

with positive energy solutions as in Eq. |T4| , m* = m + Lfs_,, and 

(pO)*=p°-U°,(p). (30) 

The subscript "i" signifies that these parameters are different for protons and neutrons. 
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As in the symmetric matter case jBrockmann & M achleidt"1990'l, evaluating the expectation 
value of Eq. ( [27| leads to a parametrization of the single particle potential for protons and 
neutrons (Eqs.i 2"4|25 ) in terms of the constants Us^ and Uy ■ which is given by 

UM = ^< PHiP)\P >= -^Us, + U°y,. (31) 

Also, 

Ur{v)= E E G,j{p,p'), (32) 

which, along with Eq. pT) , allows the self-consistent determination of the single-particle 
potential as explained below. 

The kinetic contribution to the single-particle energy is 

, , m* , , m.-m^ -|- rp- 

T,ip) = ^ <P\7-p + m\P>= ' , (33) 

and the single-particle energy is 

e*{p) = Tiip) + U,{p) = E* + I/O . (34) 

The constants m* and 

Uo,, = Us, + LfO (35) 

are convenient to work with as they facilitate the connection with the usual non-relativistic 
framework l |Haften & Tabakin||T970t . 

Starting from some initial values of m* and Uq i, the G-matrix equation is solved and a first 
approximation for Lf, (p) is obtained by integrating the G-matrix over the appropriate Fermi 
sea, see Eq. |32} . This solution is again parametrized in terms of a new set of constants, 
determined by fitting the parametrized (J,-, Eq. |3T| , to its values calculated at two momenta, 
a procedure known as the "reference spectrum approximation". The iterative procedure is 
repeated until satisfactory convergence is reached. 

Finally, the energy per neutron or proton in nuclear matter is calculated from the average 
values of the kinetic and potential energies as 

Si = j<Ti> +^ < U, > -m . (36) 
The EoS, or energy per nucleon as a function of density, is then written as 

e{p„,pp) = - (37) 

or 

e{kF,a) = f- . (38) 

Clearly, symmetric nuclear matter is obtained as a by-product of the calculation described 
above by setting a=0, whereas a=l corresponds to pure neutron matter. 
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3.2.2 Microscopic predictions of the EoS and related quantities 

In Fig. |3] we show EoS predictions for symmetric matter (solid red) and neutron matter 
(dashed black) as obtained from the Idaho calculation described in the previous section. The 
EoS from DBHF can be characterized as being moderately "soft" at low to medium density and 
fairly "stiff" at high densities. The predicted saturation density and energy for the symmetric 
matter EoS in Fig.jsjare equal to 0.185 fm^^ and -16.14 MeV, respectively, and the compression 
modulus is 252 MeV. 

The increased stiffness featured by the DBHF EoS at the higher densities originates from 
the strongly density-dependent repulsion characteristic of the Dirac-Brueckner-Hartee-Fock 
method. In Ref. (Klahn et al.[ [2006) , it is pointed out that constraints from neutron star 
phenomenology together with flow data from heavy-ion reactions suggest that such EoS 
behavior may be desirable. 

The pressure as a function of density, as discussed in the next section, plays the crucial role in 
building the structure of a neutron star. In Fig.|4]we show the pressure in symmetric matter 
as predicted by the Idaho calculation compared with constraints obtained from flow data 
jDanielewicz et al^|2002^ . The predictions are seen to fall just on the high side of the constraints 
and grow rather steep at high density. 

We show in Fig.|5]the pressure in neutron matter (red curve) and /5-equilibrated matter (green) 
as predicted by DBHF calculations. The pressure contour is again from Ref. |Danielewicz et 
[aLl [2002) . 

Next we move on to the symmetry energy as defined from Eq. |6|. In Fig.|6] we display the 
Idaho DBHF prediction for the symmetry energy by the solid red curve. The curve is seen 
to grow at a lesser rate with increasing density, an indication that, at large density, repulsion 
in the symmetric matter EoS increases more rapidly relative to the neutron matter EoS. This 
can be understood in terms of increased repulsion in isospin zero partial waves (absent from 
neutron matter) as a function of density. Our predicted value for the symmetry pressure L 
(see Eq. |8f , is close to 70 MeV. 

The various black dashed curves in Fig.|6]are obtained with the simple parametrization 

esym = C{p/pQy , (39) 

with 7 increasing from 0.7 to 1.0 in steps of 0.1, and C w 32 MeV. Considering that all of the 
dashed curves are commonly used parametrizations suggested by heavy-ion data jLi & Chen 
2005) , Fig. |6] clearly reflects our limited knowledge of the symmetry energy, particularly, but 
not exclusively, at the larger densities. 

As already mentioned in Section 2, from the experimental side intense effort is going on 
to obtain reliable empirical information on the less known aspects of the EoS. Heavy-ion 
reactions are a popular way to seek constraints on the symmetry energy, through analyses 
of observables that are sensitive to the pressure gradient between nuclear and neutron matter. 

We close this Section with demonstrating the approximately linear dependence on the 
asymmetry parameter of the single-nucleon potentials in lANM as anticipated in Eq. (To) . 
We recall that this isospin splitting is the crucial mechanism that separates proton and 
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0.5 1 1.5 

Fig. 3. DBHF predictions for the EoS of symmetric matter (solid red) and neutron matter (dashed black). 
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Fig. 4. Pressure in symmetric matter from the Idaho DBHF calculation. The shaded area corresponds to 
the region of pressure consistent with the flow data analysed in Ref. jPanielewicz et alT[|2002| . 

neutron dynamics in lANM. In Fig. [7| we display predictions obtained with three different 
NN potentials based on the one-boson-exchange model, Bonn A, B, and C ( Machleidt 1989 1. 
These three models differ mainly in the strength of the tensor force, which is mostly carried 
by partial waves with isospin equal to (absent from pure neutron matter) and thus should 
fade away in the single-neutron potential as the neutron fraction increases. In fact, the figure 
demonstrates reduced differences among the values of Lf„ predicted with the three potentials 
at large ol. 
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P (fm-3) 



Fig. 5. Pressure in neutron (red curve) and baryon-lepton (green curve) matter from the Idaho DBHF 
calculation. The shaded area corresponds to the region of pressure consistent with flow data and the 
inclusion of strong density dependence in the asymmetry term JDanielewicz et al.||2002J. 




Fig. 6. DBHF prediction for the symmetry energy (solid red) compared with various phenomenological 
parametrizations (dashed black), as explained in the text. 
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0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 



Fig. 7. The neutron and proton single-particle potentials as a function of the asymmetry parameter at fixed 
average density and momentum equal to the average Fermi momentum, which is chosen to be 1.4 fm^^. 
The red, green, and blue lines represent the predictions from the Bonn A, B, and C potentials, respectively. 

4. Neutron stars 

4.1 A brief review of basic structure equations and available constraints 

Fusion reactions in stars give rise to elements, but at the same time exhaust the nuclear fuel. 
After the fuel is exhausted, the star can die through four possible channels: it can become 
a black hole, a white dwarf, a neutron star, or it can disassemble completely. The ultimate 
outcome depends on the mass of the original star. If the mass is larger than about four solar 
masses, the star may become a supernova which, in turn, may result either in a neutron star 
or a black hole. 

Neutron stars contain the most dense form of matter found in the universe and, therefore, are 
unique laboratories to study the properties of highly compressed (cold) matter. They are held 
together by gravity and neutron degeneracy pressure. (In contrast, white dwarfs are kept in 
hydrostatic equilibrium by gravity and electron degeneracy pressure.) 

Although neutron stars were predicted as early as in the 1930's, hope for their observation 
remained slim for a long time. In 1967, strange new objects, outside the solar system, were 
observed at the University of Cambridge. They were named pulsars, as they emitted periodic 
radio signals. To date, about 1700 pulsars have been detected, many in binary systems. 

Typically, detection of thermal radiation from the surface of a star is the way to access 
information about its properties. Furthermore, the possibility of exploring the structure of 
neutron stars via gravitational waves makes these exotic objects even more interesting. 
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The densities found in neutron stars range from the density of iron to several times normal 
nuclear density. Most of the mass consists of highly compressed matter at nuclear and 
supernuclear densities. The surface region is composed of normal nuclei and non-relativistic 
electrons, with typical mass densities in the range 10* < e < 10^ g cm^-^. As density increases, 
charge neutrality requires matter to become more neutron rich. In this density range (about 
g cm neutron-rich nuclei appear, mostly light metals, while electrons 
become relativistic. This is the outer crust. Above densities of approximately 10^^ g cm^^, free 
neutrons begin to form a continuum of states. The inner crust is a compressed solid with a 
fluid of neutrons which drip out and populate free states outside the nuclei, since those have 
become neutron-saturated. Densities in the inner crust range between 10^^ and 10^* g cm^''. 
At densities equal to approximately 1 /2 of saturation density, clusters begin to merge into a 
continuum. In this phase, matter is a uniform fluid of neutrons, protons, and leptons. Above 
a few times nuclear matter density, the actual composition of stellar matter is not known. 
Strange baryons can appear when the nucleon chemical potential is of the order of their rest 
mass. Meson production can also take place. At even higher densities, transitions to other 
phases are speculated, such as a deconfined, rather than hadronic, phase. The critical density 
for such transition cannot be predicted reliably because it lies in a range where QCD is non 
perturbative i Se drakian[|2007^ . 

The possibility has been speculated that the most stable state at zero pressure may be ii, d, s 
quark matter instead of iron. This would imply that strange quark matter is the most stable 
(in fact, the absolutely stable) state of strongly interacting matter, as originally proposed by 
Bodmer jBodmer||1971| , Witten l |Witten[[l984) , and Terazawa jTerazawa l. 1989). In such case, 
hyperonic and hybrid stars would have to be metastable with respect to stars composed 
of stable three-flavor strange quark matter (Weber, 1999), which is lower in energy than 
two-flavor quark matter due to the extra Fermi levels open to strange quarks. Whether or 
not strange quark stars can give rise to pulsar glitches (which are observed sudden small 
changes in the rotational frequency of a pulsar), may be a decisive test of the strange quark 
matter hypothesis ( |Weber[|1999| . 

The maximum gravitational mass of the star and the corresponding radius are the typical 
observables used to constraint the EoS. The gravitational mass is inferred mostly from 
observations of X-ray binaries or binary pulsars. Determination of the mass provide a unique 
test of both theories of nuclear matter and general relativity. The pulsar in the Hulse-Taylor 
binary system has a mass of 1.4408± 0.0003 Mr:,, to date the best mass determination. 

At this time, one of the heaviest neutron stars (with accurately known mass) has a mass of 
1.671 ± 0.008 Mq l |Champion et al. 2008^ . The observation of an even heavier star has been 



confirmed recently, namely J1614-2230, with a mass of 1.97±0.04 M© pemorest et al.[ |2010) 
This value is the highest yet measured with this certainty and represents a challenge for the 
softest EoS. We also recall that an initial observation of a neutron star-white dwarf binary 
system had suggested a neutron star mass (PSR J0751+1807) of 2.1±0.2M:3 ( Nice et al.[[2005^ . 
Such observation, which would imply a considerable constraint on the high-density behavior 
of the EoS, was not confirmed. 



The minumum mass of a neutron star is also a parameter of interest. For a cold, stable system, 
the minimum mass is estimated to be 0.09 M© l |Lattimer & Prakash| [2007) . The smallest 
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reliably estimated neutron star mass is the companion of the binary pulsar J1756-2251, which 
has a mass of 1.18± 0.02 M© ^Faulkner et"aL][2004^ . 



Measurements of the radius are considerably less precise than mass measurements ^Lattimer 



& Prakash) | 2007| . No direct measurements of the radius exist. Instead, the observed X-ray 



flux, together with theoretical assumptions (Weber 1999), can provide information on the 
radiation or photospheric radius, Roa, which is related to the actual stellar radius by Rco = 
R{1 — 2GM/Rc-^)^^^^. Estimates are usually based on thermal emission of cooling stars, 
including redshifts, and the properties of sources with bursts or thermonuclear explosions 
at the surface. A major problem associated with the determination of radii is that the distance 
from the source is not well known, hence the need for additional assumptions. Much more 
stringent constraints could be imposed on the EoS if mass and radius could be determined 
independently from each other. 

Another bulk property of neutron stars is the moment of inertia, I. For softer EoS, both mass 
and radius are smaller and so is I. From observations of the Crab nebula luminosity, a lower 
bound on the moment of inertia was inferred to be I >4-8 x 10** g cm^, see Ref. ^ Weber,, 1999| 
and references therein. A measurement of the moment of inertia within 10%, together with 
the information on the mass, would be able to discriminate among various EoS ('Lattimer & 



[Prakas h, 2007 1. To date, the best determination of the moment of inertia is the one for the Crab 



pulsar |Beiger & Haensel 2003^ which would rule out only very soft EoS i Lattimer & Prakash 
[2007) . 

A proton-neutron star is the result of a supernova explosion resulting from the gravitational 
collapse of a massive star core. Nearly all of the remaining binding energy is carried away 
by neutrinos during the first few tens of seconds of the evolution. Thus neutrino emission is 
very efficient as a cooling mechanism, with the internal temperature dropping to about 10^" K 
within a few days. Cooling through neutrino emission continues for a long time (in the order 
of 1,000 years), until the temperature drops to about 10* K, at which point photon emission 
becomes the dominant cooling mechanism. Neutrino luminosity and emission timescale are 
controlled by several factors including the total mass of the (proton-neutron) star and the 
opacity of neutrinos at high densities, which is sensitive to the EoS of dense hadronic matter. 

Gravitational waves are a less conventional way to probe neutron star properties. Compact 
stars in binary systems are epected to produce gravitational radiation. In turn, emission of 
gravitational waves causes decay of the mutual orbits and eventually merger of the binary 
system. Because of the merger timescale (250 million years for PSR B1913+16, for instance, 
and 85 million years for PSR J0737-3039), it can be expected that many such decaying binary 
systems exist in the galaxy and emit large amounts of gravitational radiation. The observation 
of gravitational waves has the potential to set strong constraints on masses and radii (see 
{Lattimer & Prakash[[2007) and references therein). 

A theoretical estimate of the maximum possible mass of a neutron star was performed by 
Rhoades and Ruffini ( Rhoades & Ruffim} 1974) on the following assumptions: 1) General 



relativity is the correct theory of gravitation; 2) the EoS satisfy the Le Chatelier's principle 
{dP/de > 0) and the causality condition, dP/de < c^; and 3) the EoS below some matching 
density is known. On this basis, they determined that the maximum mass of a neutron 
star cannot exceed 3.2 solar masses. Abandoning the causality condition, which would hold 
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exactly only if stellar matter is neither dispersive nor absorptive, this limit can be as high as 5 
solar masses due to the increased stiffness of the EoS at supernuclear densities. 

The maximum mass and the radius of a neutron star are sensitive to different aspects of the 
EoS. The maximum mass is mostly determined by the stiffness of the EoS at densities greater 
than a few times saturation density. The star radius is mainly sensitive to the slope of the 
symmetry energy. In particular, it is closely connected to the internal pressure (that is, the 
energy gradient) of matter at densities between about l.dpo and 2-3po jLattimer & Prakash 



2007| . Non-nucleonic degrees of freedom, which typically make their appearance at those 



densities, are known to have a considerable impact on the maximum mass of the star. The 
latter is predicted by the equation of hydrostatic equilibrium for a perfect fluid. 

In general relativity, the invariant interval between two infinitesimally close space-time events 
is given by 

ds^ = g^lidx^dx^ , (40) 

where g^^ is the space-time metric. For a spherically symmetric space-time, the most general 
static line element consistent with all required synnmetries has the form 

ds^ = -f{r)dt^ + gir)dr^ + h^{r)(sin^e dip^ + dfi) . (41) 

Choosing a radial coordinate r such that = yields the so-called "standard" form of the 
metric. 

The equation of hydrostatic equilibrium (the TOV equation) determines the form of the metric 
functions along with the pressure and the total mass-energy density as a fimction of the radial 
coordinate in the interior of the star. It reads 

dP{r) _ _G {P{r)+e{r)){M{r)+4:TTpP{r)/c^) 

'liT ~ (? r(r-2GM(r)/c2) ' ^ ' 

with 

^=47rr2e(r), (43) 
where e is the total mass-energy density. The star gravitational mass is 

M(R) = r Anr^eMdr , (44) 
Jo 

where R is the value of r where the pressure vanishes. It's worth recalling that no mass limit 
exists in Newtonian gravitation. 

The pressure is related to the energy per particle through 

P(p)=p^^. (45) 

The structure equations of rotationally deformed compact stars are much more complex 
than those of spherically symmet ric stars jWeber| 1999} presented here. The most rapidly 



rotating pulsar, PSR J1748-2446 (Hesse ls et al.| |2006), isbelieved to rotate at a rate of 716 Hz, 
although an X-ray burst oscillation at a frequency of 1122 Hz was reported jKaaret et al.[|2006"| , 
which may be due to the spin rate of a neutron star. Naturally, the maximum mass and the 
(equatorial) radius become larger with increasing rotational frequency. 
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4.2 Composition of /3-stable matter 

Assuming that only neutrons, protons, and leptons are present, the proton fraction in stellar 
matter under conditions of /3-equilibrium is calculated by imposing energy conservation 
and charge neutrality. The resulting algebraic equations can be foimd in standard literature 
jGlendenningl [T997| . The contribution to the energy density from the electrons is written as 

= ^{3n^p,)^/\ (46) 

whereas for muons we write 

e,, = p,,m,,c2 + {hc)^ \ . (47) 

These contributions are added to the baryonic part to give the total energy density. The 
derivative of the total energy per particle with respect to the fraction of a particular species is 
the chemical potential of that species. The conditions 

fip + }h = f'n ; = Ve ; Pp = Pfi + Pe , (48) 

allow to solve for the densities (or fractions) of protons, electrons, and muons. Near the 
saturation density, when the muon fraction is close to zero, one can estimate the equilibrium 



proton fraction, Xp, to be ' Lattimer & Prakash 



2007 



/4£'sy,„(|Oo)N|3 



/{3n^Po)- (49) 



The fractions of protons, electrons, and muons as predicted with the DBHF equation of state 
are shown in Fig. |8] The critical density for the proton fraction to exceed approximately 1/9 
and, thus, allow cooling through the direct Urea processes, 

n ^ p + e + Ve and p + e — >■ « + , (50) 

is about 0.36 — 0.39 fm^^. Notice that, due to the relation between symmetry energy and 
proton fraction, large values of the symmetry energy would make the star cool rapidly. In 
fact, already in earlier studies ( BogutaJ |198 1| tha rapid cooling of neutron stars and the 
corresponding high neutrino luminosity was understood in terms of neutron fi decay and 
large proton fractions. 

At densities close to normal nuclear density, protons and neutrons are the only baryonic 
degrees of freedom. As density increases, other baryons begin to appear, such as strange 
baryons or isospin 3/2 nucleon resonances. Hyperonic states can be classified according to 
the irreducible representation of the S!i(3) group. The octet of baryons that can appear in 
neutron matter includes nucleons. A, YP'^, and 'SP'^ . 

Neglecting the nucleon-hyperon interaction, the threshold for stable hyperons to exist in 
matter is determined by comparing the hyperon mass with the neutron Fermi energy, which 
is the largest available energy scale in neutron-rich matter. We consider cold neutron stars, 
after neutrinos have escaped. Strange baryons appear at about 2-3 times normal density 
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jBaldo et"aL 1998| , an estimate which is essentially model independent, through the processes 
M + w — >■ p + Yr and n + w — )• n + A. The equilibrium conditions for these reactions are 



2}in = }lp + }ll.- ; fin = fiA- 

Also, we have 

the equations above being special cases of 

fi = bfi„- qjie, 



(51) 
(52) 

(53) 



where b and q are the baryon number and the charge (in units of the electron charge) of the 
particular species with chemical potential pi. Together with the charge neutrality condition 
and baryon number conservation. 



Pp = pe + p,i + PL- ; P = Pii + Pp + Pi.- + Pa , 
the above system allows to determine the various particle fractions. 



(54) 



Naturally, the composition of matter at supra-nuclear densities determines the behavior of 
stellar matter. It is also speculated that a transition to a quark phase may take place at very 
high densities, the occurrence of which depends sensitively on the properties of the EoS in 
the hadronic (confined) phase. The presence of hyperons in the interior of neutron stars is 
reported to soften the equation of state, with the consequence that the predicted neutron star 
maximum masses become considerably smaller jSchulze et al. 2006|. Strange baryons are not 
included in the predictions shown below. 
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Fig. 8. Proton (solid black), electron (dashed red), and muon (dotted blue) fractions in /J-stable matter as 
a function of total baryon density as predicted by the DBHF model. 
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Fig. 9. The baryon density and the mass-energy density profile for a neutron star with the maximum mass 
allowed by each EoS model. 

4.3 Microscopic predictions of neutron star properties 

We are now ready to move to applications of our EoS to compact stars. 

As explained in Section 3, the DBHF model does not include three-body forces explicitely, 
but effectively incorporates the class of TBF originating from the presence of nucleons and 
antinucleons (the "Z-diagrams" in Fig.|2]l, see discussion in Section 3.1. In order to broaden our 
scopes, we will compare our predictions with those of other microscopic models. As the other 
element of our comparison, we will take the EoS from the microscopic approach of Ref. jLi et 



al. 2008). There (and in previous work by the same authors), the Brueckner-Hartree-Fock 



(BHF) formalism is employed along with microscopic three-body forces. In particular, 
in Ref.l jLi & Schuize} [2008^ the meson-exchange TBF are constructed applying the same 
parameters as used in the corresponding nucleon-nucleon potentials, which are: Argorme 
V18 jWiringa et aI.|[T995 ) (V18), Bonn B lMachl eidt||T989) (BOB), Nijmegen 93 (St ocks et af 



1994) (N93). The popular (but phenomenological) Urbana TBF ^Fieper et al., ,2001| (UIX) is 



Will-be-set-by-IN-TECH 



21 





Fig. 10. The baryon density and the mass-energy density profile for a neutron star with a mass of 1.4 solar 
masses. 

also utilized in Ref. jLi & Schulze[[2008) . Convenient parametrizations in terms of simple 
analytic functions are given in all cases and we will use those to generate the various EoS. We 
will refer to this approach, generally, as "BHF + TBF". 

At subnuclear densities all the EoS considered here are joined with the crustal equations of 
state from Harrison and Wheeler ((Harrison et al. 19651 and Negele and Vautherin iNegele &] 
|Vautherirt}|1973^. The composi tion of the crust is crystalline, with light ^Harrison et al.,|1965] l 
or heavy jNegele & Vautherin |1973^ metals and electron gas. 



We begin with showing the baryon number density and the mass-energy density profile of the 
star, see Fig.|9] For each EoS model, the maximum mass configuration is considered. Thus the 
models differ in their central density, which, in turn, impact the radius. The relations shown 
in Fig.|9]are insightful, as they reveal the detailed structure of the star at each radial position. 
Furthermore, the compactness of the star, whose density profile is reminescent of the one in a 
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Fig. 11. Pressure profile for the maximum-mass star allowed by each EoS model. 
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Fig. 12. Pressure profile predicted by the various models for a 1.4 solar mass star 

nucleus, a system 55 orders of magnitude lighter, is apparent. In Fig.jlO] the same quantities 
are shown for a star with a mass of 1.4 solar masses, the most probable mass of a neutron star. 

The models labeled as UIX and BOB have the smallest and largest radius, respectively, as 
can be seen from the figure. We also see that the star's outer regions, that is, for energy 
densities less than about lO^^g cm^'', are influenced quite strongly by differences in the 
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Fig. 13. Energy density us. baryon number density for the various EoS being addressed in the text. The 
maximum-mass model is considered in each case. 

various EoS models. Note that the UIX model, with the smallest radius, can tolerate larger 
central densities. 

Of interest is also the pressure profile for the maximum-mass star in each model, which is 
shown in Fig.llllfor the maximum mass and in Fig. 12 for a 1.4 solar mass star. 



The e(r) and p{r) relations from Fig.[9]are combined to provide the £{p) relation within the 
star as shown in Fig. 13 for the maximum mass. Again, we see that the stiffest (BOB) and 
softest EoS (UIX) support the smallest and largest central densities, respectively. At the same 
time, these two EoS predict the largest (BOB) and smallest (UIX) maximum mass, see below. 



In Fig. 14 we show the mass-radius relation for a sequence of static neutron stars as predicted 
by the various models. All models besides DBHF share the same many-body approach 
(BHF+TBF) but differ in the two-body potential and TBF employed. The differences resulting 
from the use of different NN potentials can be larger than those originating from emplying 
different many-body approaches. This can be seen by comparing the DBHF and BOB curves, 
both employing the Bonn B interaction (although in the latter case the non-reltivistic, r-space 
version of the potential is adopted). Overall, the maximum masses range from I.SMq (UIX) to 
2.5M0 (BOB). Radii are less sensitive to the EoS and range between 10 and 12 km for all models 
under consideration, DBHF or BHF+TBF. Concerning consistency with present constraints, 
the observations reported in Section 4.1 would appear to invalidate only the model with 
the smallest maximum mass, UIX. Notice, further, that phenomena such as condensation of 
mesons may soften the EoS considerably at supernuclear density as condensation would bring 
loss of pressure. 
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Also of interest is the star baryon number, A, which is obtained by integrating the baryon 
density over the proper volume jWeber 1999}. Namely, 



A = An r dr , ^''^^ . (55) 

■/o v'l-2GM(r)/(rc2) 

Defining the star's baryon mass as 

Ma = m,! A , (56) 
where m„ is the mass of the baryon, one calculates the star's binding energy as 

Eb = M-Ma. (57) 

The baryon number and the star binding energy as a function of the central density (in units 
of nuclear matter density) are shown in Fig. [15] and Fig. [16] for the various models. We see 
that the baryon number for stable stars is approximately equal to 10^^ — 10^''. A much higher 
value would make the star unstable with respect to gravitational collapse. The binding energy 
displayed in Fig. 16 is defined in units of the solar mass. From the definition, Eq.|57] negative 
values signify a bound system. Typically, the binding energy changes the sign for masses less 
than 0.1 solar masses. The binding energy is a potentially observable quantity, since neutrinos 
from a supernova carry information about the gravitational binding energy of the neutron star 
that has resulted from the explosion. 

Next, we calculate the gravitational redshift predicted by each model. The redshift is defined 
as 

z = ^ - 1 , (58) 

where Ve and Voo are the photon frequencies at the emitter and at the infinitely far receiver. The 
photon frequency at the emitter is the inverse of the proper time between two wave crests in 
the frame of the emitter, 

^ = i-gf^ix'^dx'')-^^^ , (59) 
with a similar expression for the frequency at the receiver. Then 

ii-gooV^'dx^U ■ ^ ' 

Assuming a static gravitational field, in which case the time dx'^ between two crests is the 
same at the star's surface and at the receiver, and writing ggg as the metric tensor component 
at the surface of a nonrotating star yield the simple equation 

Notice that simultaneous measurements of Rca and z determines both R and M, since 

R = Roo{l + z)-\ (62) 

and ^ 

M= ^R^{l + z)-\l-{l + z)-^]. (63) 
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Fig. 14. Static neutron star mass-radius relation for the models considered in the text. 

In Fig.[l7|we show the gravitational redshift as a function the mass for each model. Naturally 
the rotation of the star modifies the metric, and in that case different considerations need to 
be applied which result in a frequency dependence of the redshift. We will not consider the 
general case here. 

We conclude this section with showing a few predictions for the case of rapidly rotating stars. 
The model dependence of the mass-radius relation is shown in Fig.[T8| The 71 6 Hz frequency 
corresponds to the most rapidly rotating pulsar, PSR J1748-2446, jHessels et al.|[2006| although 
recently an X-ray burst oscillation at a frequency of 1122 Hz has been reportedfKaaret et al.| 
[2006) which may be due to the spin rate of a neutron star. As expected, the maximum mass 
and the (equatorial) radius become larger with increasing rotational frequency. 

In Fig. [19] we show the moment of inertia at different rotational speeds (again, for all models). 
These values are not in contradiction with observations of the Crab nebula luminosity, from 
which a lower bound on the moment of inertia was inferred to be I >4-8 x 10'** g cm^, see 
Ref. jWeberl|1999) and references therein. 

Clearly, at the densities probed by neutron stars the model dependence is large, but presently 
available constraints are still insufficient to discriminate among these EoS. The model 
dependence we observe comes from two sources, the two-body potential and the many-body 
approach, specifically the presence of explicit TBF or Dirac effects. The dependence on the 
two-body potential is very large. Typically, the main source of model dependence among NN 
potentials is found in the strength of the tensor force. Of course, differences at the two-body 
level impact the TBF as well, whether they are microscopic or phenomenological. 

5. Polarized isospin-asymmetric matter 

Before concluding this chapter, we like to touch upon the issue of polarization in lANM. 
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Fig. 15. Baryon number as a function of the central density (in units of eg = 2.5 10^* g cm ^) for tfie 
models considered in the text. 
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Fig. 16. Binding energy (in units of the solar mass) as a function of the central density for the EoS models 
considered in the text. 

When both isospin and spin asymmetries are present, constraints are much more difficult to 
obtain and predictions regarding magnetic properties of nuclear matter are sometimes found 
to be in qualitative disagreement with one another. This is especially the case with regard to 
the possibility of spontaneous phase transitions into spin ordered states, ferromagnetic (FM, 
with neutron and proton spins alligned), or antiferromagnetic state (AFM, with opposite spins 
for neutrons and protons). Notice that the presence of polarization would impact neutrino 
cross section and luminosity, resulting into a very different scenario for neutron star cooling, 
which is why we find it appropriate to briefly discuss this issue here. 



Recently, we have extended the framework described in Section 3 > Sammarruca 
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Will-be-set-by-IN-TECH 




12 3 



Mass (M,„„) 



Fig. 17. Gravitational redshift for all models. For each model, the corresponding sequence of static stars 
is considered. 

include both spin and isospin asymmetries of nuclear matter and calculated the energy per 
particle imder extreme conditions of polarization. The existence (or not) of a possible phase 
transition can be argued by comparing the energies of the fully polarized and the impolarized 
phases. 

In a spin-polarized and isospin asymmetric system with fixed total density, p, the partial 
densities of each species are 

pH = Pnu + Pnd ■ Pp = Ppu +Ppd , P = Pn + Pp , (64) 

where u and d refer to up and down spin-polarizations, respectively, of protons (p) or neutrons 
(n). The isospin and spin asymmetries, a, /3„, and /3p, are defined in a natural way: 

Pn - pp Pnu - Pnd n PP" ~ Ppd 
OL = — — , /3„ = , /3p = . (65) 

P Pn Pp 



The single-particle potential of a nucleon in a particular tu state, (Jto-, is now the solution of a 
set of four coupled equations, which are the appropriate extension of Eqs. j24|25^. They read 





— ^nu,nu 




H" ^nii,pii 


~l~ ^uu,pd 


(66) 




^nd,nu 


~^ ^nd,nd 


~l~ ^nd,pu 


+ ^nd,pd 


(67) 




— Upu,nu 


~^ ^pu,nd 


H~ Upu,pu 


~l~ ^pii,pd 


(68) 



^pd ^pd,nu ^pd,nd ~^ ^pd,pu ~t~ ^pd,pd / (69) 

to be solved self-consistently along with the two-nucleon G-matrix. In the above equations, 
each Ujcr,T'cr' term contains the appropriate (spin and isospin dependent) part of the 
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Fig. 18. Mass-radius relation for the models considered in the text and for different rotational frequencies, 
interaction, Grcr/T'cr'- More specifically, 

Ura{k)= ELL <rcT,T'cT'\G{k,q)\Ta,T'a' >, (70) 

(r'=u,d T'=n,p cjKk}''^' 
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Fig. 19. Moment of inertia for the models considered in tfie text and for different rotational frequencies, 
where the third summation indicates integration over the Fermi seas of protons and neutrons 
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Fig. 20. The energy per particle as a function of density and variuos degrees of proton and neutron 
polarizations in symmetric matter (left) and asymmetric matter (right). In both frames, the (blue) 
dashed line corresponds to totally polarized neutrons and unpolarized protons (/5«=1, fip=0); the (green) 
dash-dotted line is the prediction for the FM state (/3h=1, /3p=l); the (red) dotted line shows the energy of 
the AFM state (/in=l/ /5p=-l)- The (black) solid line shows the predictions for unpolarized matter 



with spin-up and spin-down, 
presence of spin polarization. 



Notice that this equation is the extension of Eq. |32| in the 



In the left panel of Fig.|20] we show, in comparison with unpolarized symmetric matter (solid 
line): the EoS for the case of fully polarized neutrons and completely unpolarized protons 
(dashed line); the EoS for the case of protons and neutrons totally polarized in the same 
direction, that is, matter in the ferromagnetic (FM) state ( dashed-dotted line); the EoS for 
the case of protons and neutrons totally polarized in opposite directions, namely matter in the 
antiferromagnetic (AFM) state ( dotted line). A similar comparison is shown in the right panel 
of Fig. l |20| l, but for isospin asymmetric matter. (Notice that all predictions are invariant imder 
a global spin flip.) 

We conclude that, for both symmetric and asymmetric matter, the energies of the FM and 
AFM states are higher than those of the corresponding unpolarized cases, with the AFM state 
being the most energetic. Thus, a phase transition to a spin-ordered state is not anticipated 
in our model. This conclusion seems to be shared by predictions of microscopic models, such 
as those based on conventional Brueckner-Hartree-Fock theory ( Vidafia & Bombaci} |2002V 
On the other hand, calculations based on various parametrizations of Skyrme forces result 
in different conclusions. For instance, with the SLy4 and SLy5 forces and the Fermi liquid 
formalism a phase transition to the AFM state is predicted in asymmetric matter at a critical 
density equal to about 2-3 times normal density (Ilsayev & Yang |2004| . 



In closing this brief section, it is interesting to remark that models based on realistic 
nucleon-nucleon potentials, whether relativistic or non-relativistic, are at least in qualitative 
agreement with one another in predicting more energy for totally polarized states (FM or 
AFM) up to densities well above normal density. 
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6. Summary and Conclusions 

In this chapter, we have been concerned with the nuclear equation of state of isospin 
asymmetric nuclear matter, the main iaput for calculations of the properties of compact stars 
as well as a variety of other systems, such as the neutron skin of neutron-rich nuclei. 

After describing our microscopic approach to the development of the equation of state for 
nuclear matter and neutron-rich matter, we presented a brief review of the structure equations 
leading to the prediction of neutron star properties. Microscopic predictions from different 
models employing three-body forces along with the non-relativistic Brueckner-Hartree-Fock 
method have also been shown for comparison. Large model dependence is seen among 
predictions, especially those involving the highest densities. 

Rich and diverse effort is presently going on to improve the available constraints on the EoS 
or find new ones. These constraints are usually extracted through the analysis of selected 
heavy-ion collision observables. At the same time, partnership between nuclear physics and 
astrophysics is becoming increasingly important towards advancing our understanding of 
exotic matter. The recently approved Facility for Rare Isotope Beams (FRIB), thanks to new 
powerful technical capabilities, will forge tighter links between the two disciplines, as it will 
allow access to rare isotopes which play a critical role in astrophysical processes but have not 
yet been observed in terrestrial laboratories. 
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